Gutzwiller theory of band magnetism in LaOFeAs 
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We use the Gutzwiller variational theory to calculate the ground-state phase diagram and quasi- 
particle bands of LaOFeAs. The Fe3d-As4p Wannier-orbital basis obtained from density-functional 
theory defines the band part of our eight-band Hubbard model. The full atomic interaction between 
the electrons in the iron orbitals is parameterized by the Hubbard interaction U and an average 
Hund's-rule interaction J. We reproduce the experimentally observed small ordered magnetic mo- 
ment over a large region of {U, J) parameter space. The magnetically ordered phase is a stripe 
spin-density wave of quasi-particles. 

PACS numbers: 71.10.Fd, 71.20.Be, 71.27.+a 



LaOFeAs is the first high-temperature Fe-based super- 
conductor [l[ . Its crystal structure is made of alternating 
layers of (LaO)+ and (FeAs)". Below = 150 K, the 
material undergoes a structural phase transition from 
tetragonal to orthorhombic, follovifcd by the onset of a 
stripe spin-density wave (SDW) order in the FeAs layers 
with wave vector q = (0, tt). The values of the ordered 
moments vary from m — 0.3^b to 0.8^b 0-11] • Under 
pressure or upon doping, the electrons in the FeAs layers 
become superconducting at Tc = 56 K [1, 0] . The pres- 
ence of Fe-pnictogen or chalcogen layers and the prox- 
imity of magnetism and superconductivity with fairly 
high Tc are common to almost all the iron-based super- 
conductors known to date, and place these compounds 
in the class of unconventional superconductors like the 
cuprates. 

However, the cuprates and pnictides show important 
differences. In contrast to the cuprates, the normal 
phase of the pnictides is always metallic; in fact, density- 
functional theory (DFT) calculations provide a qualita- 
tively correct description of the electronic structure and 
Fermi surface (FS) of the paramagnetic (PM) phase 0,[1|, 
and reproduce the correct magnetic structure of the 
SDW state [9]. However, sizable renormalizations of the 
DFT quasi-particle band masses are observed in optics, 
angular- resolved photo-emission spectroscopy (ARPES), 
and de-Haas-van- Alphen measurements [5[. Moreover, 
the values of the ordered moments calculated by DFT 
in the SDW phase, todft ~ 1-8/iB or higher, are much 
larger than the experimentally observed values 

This discrepancy has lead to an intense ongoing debate 
on the nature of magnetism (localized vs. itinerant), and 
on the mechanism of suppression of the magnetic mo- 
ment (long-range spin fluctuations vs. strong local elec- 
tronic correlations). Adopting the correlated-electron 
viewpoint l3, Hi, one faces the problem that very few 
theoretical methods are available for the description of 
itinerant magnetism in multi-band systems. 



In this work, we employ the Gutzwiller variational the- 
ory (GT) which describes the ground state and quasi- 
particle excitations of multi-band Fermi liquids to study 
the magnetic phase diagram of LaOFeAs; GT treats 
the full atomic interactions and is numerically much 
less involved than, e.g.. Dynamical Mean-Field Theory 
(DMFT) [iMl- This enables us to resolve the smah 
energy differences between the PM and SDW phases in 
the pnictides. We shall find that the striped SDW in 
LaOFeAs has a small ordered moment over a large re- 
gion of (C/, J) parameter space. This small-moment phase 
can be understood as the band magnetism of correlated 
quasi-particles. 

We investigate the two-dimensional eight-band Hub- 
bard model 



H = Ho + Hc 



E 



(1) 



where Hq describes the on-site energies and hopping in- 
tegrals of electrons with spin a \. between FeSd-Asdp 
Wannier orbitals h,b' on sites i, j (NMTO) [8]. The cor- 
responding DFT bandstructure is shown in Fig. [1] We 
have neglected the interlayer hoppings and have used the 
glide plane to unfold to the Brillouin zone which contains 
only one formula unit per cell [83 ■ 

The local Hamiltonian Hc,i describes the Coulomb in- 
teraction between the d-electrons. We use the Hubbard 
interaction U and an averaged Hund's-rule interaction J 
to derive the Racah coefficients A,B,C with the atomic 
ratio C/B = 4 which parameterize our multi-band inter- 
action 13|. In the following, until in Figs. S] and [5l we 
give the results for a characteristic set of values U = 8eV 
and J = 0.6 eV [{U, J) = (8, 0.6) eV]. This set, for which 
J — 0.075C/, gives an ordered moment of m = 0.74^b, 
intermediate between two recent experimental values. 

We fix the average number of 3d-electrons per Fe and 
of the 4p-electrons per As at their PM DFT- values, Ud = 
7.47 and Up = 4.53, respectively jSj. The polar-covalent 
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FIG. 1. (Color online) Fe3d-As4p bandstructure of LaOFeAs 
along the high-symmetry lines of the 2D Brillouin zone con- 
taining one formula unit per cell. Left: DFT 8] (dashed) 
and GT (solid) in the PM phase for (17, J) = (8, 0.6) eV. The 
bands are lined up at the Fermi level {E = 0), energies on the 
ordinate are in eV, the DFT bandwidth is Wdft = 6.5 eV. 
Right: GT bands for {U, J) = (8, 0.6) eV in the PM (dashed) 
and SDW (solid) phases in the Brillouin zone folded along the 
line ^FY-^XM. High symmetry points and the directions x, 
y and X, Y are defined in Fig. [S] Dominant orbital and spin 
characters are indicated. Spin-f are minority and spin-j- are 
majority spin states. For strongly covalent bands, z/xy means 
more z and less xy character and vice versa [8|. 



bonding between Fe and As thus transfers 1.47 electrons 
from As to Fe compared with the ionic picture (dfip^) 
implicit in the often used five-band model [13] • It has 
been shown previously that the ionic picture does 
not lead to a satisfactory description of the magnetically 
ordered phase of LaOFeAs within GT. We do not fix the 
occupations of each of the five Sd orbitals but there is 
little charge flow between them as a function of (C/, J) 
for the parameter values considered. 

The true ground state of if in ([T]) is approximated by 
the Gutzwiller variational wave function 



(2) 



where j^'o) is a product state of Bloch orbitals, and the 
local Gutzwiller correlator is defined as 



^^ = Er^r|r)a(r| 



(3) 



Variational parameters Ar have been introduced for each 
multiplet state \T)i on the iron sites, i.e., the eigenstates 
of Hc,i- Note that |\l/o) is also a variational object which 
we determine from the minimization of the variational 
energy functional resulting from the wave functions ([2|). 

As shown in Refs. 15|, ll6j , expectation values can be 
evaluated without further approximations in the limit of 
infinite spatial dimensions. The energy functional de- 
rived in this limit is then used as an approximation for 



finite-dimensional systems. In addition, GT provides the 
Landau-Gutzwiller quasi-particle bandstructure for com- 
parison with ARPES data [13, Ell- 

We begin with the hight-temperature PM phase and 
present the quasi-particle bandstructure in Fig. [T] The 
electron-electron interaction between the d-electrons is 
seen to reduce the overall bandwidth. Filled bands with 
a large character are shifted upwards in energy and 
empty, 3d-like bands downwards. Also, the three lowest 
bands are effected since of the 6 electrons in these As p- 
like bands, 1.82 are Fe d. 





FIG. 2. (Color online) PM Fermi surface. Left: DFT; Right: 
GT for {U, J) — (8, 0.6) eV. The lower and upper case coor- 
dinate systems used to designate the orbitals are shown. 

In Fig.Hwe show the FS calculated with DFT and GT 
for ([/, J) = (8, 0.6) eV for the PM phase. In GT the cor- 
relations shift the xy band down by 130 meV and thereby 
empty the F-centered hole pocket. The concomitant low- 
ering of the Fermi level narrows the X and Y-centered 
electron pockets and makes the inner M-centered hole 
pocket expand. The lowering of the empty XY/Yz- 
like band near M steepens the lower, Fz-like conduction 
band [s] and thereby shrinks the outer M-centered hole 
pocket along M-F. The small size of these effects con- 
firms the conventional wisdom that DFT bands provide 
a reasonable description of the FS. On the other hand, 
the agreement of the DFT FS with experiment does not 
indicate that the electron-electron interactions are weak. 

The correlation-induced mass renormalization on the 
FS is defined as the renormalization of the Fermi ve- 
locity, ni* /m — wf.dft/vf. It depends on the orbital 
character of the band states involved. For pure d-bands, 
the renormalization is qd ~ 1-5 for (C/, J) = (8,0.6)eV 
and grows by 3% per eV increase of U . The mass en- 
hancement decreases towards unity with increasing p- 
character. Also band-shifts influence the effective masses 
on the Fermi surface. For the six points where the 
Fermi surface cuts the high-symmetry lines in Fig. [21 
{m*/m)n = 1.6, 1.2, 1.3, 1.3, 1.4, 1.3 for the orbitals n = 
xz/y, xyjz^ xz, yz, Xz, Yz. Note that the xz/y-like DFT 
band in F-X direction has a strong curvature in the vicin- 
ity of X and, as the correlation-induced downwards shift 
of the xy band drags the Fermi level along, it is placed 
lower in the xz/y band. This increases {m* /m),j.z/y be- 
yond the renormalization factor qd- 



3 



0.5r 

0.4 
0.3 
0.2 
0.1 
0.0 



, n(GT) 

] n(4'„) 



0123456789 10 
number of d-electrons 



0.40 
0.35 
0.30 
0.25 
0.20 
0.15 
0.10 
0.05 
0.00 



i 



1 5(GT) 



0.0 0.5 1.0 1.5 2.0 2.5 

local spin, s 



FIG. 3. (Color online) Local d-charge distribution p{n) (left) 
and local spin distribution p{s) (right) for the PM GT with 
([/, J) — (8, 0.6) eV and for the PM single-particle wave func- 
tion |*o)- 
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For ([/, J) — (8, 0.6) cV, the dominant Sd-configur- 
ations are and id^ in the optimized PM Gutzwiller 
wave function j^'^'^'). Occupancies with = 6 or = 9 
electrons are rare, and all other charge states are essen- 
tially forbidden, see Fig. [31 The absence of large charge 
fluctuations is typical for correlated electron systems. 
The probability distribution function p[s) for finding lo- 
cal spins with size < s < 5/2 in is very similar to 
the distribution in the single-particle product state {^o)- 
In GT, the correlation enhancement of the local spin mo- 
ment is small. This shows that the system is far from 
the local-moment situation where we would solely find 
atoms with Hund's-rule spins s = 3/2 for 3d^ and s = 1 
for 3d*. From the expectation values for the local spin, 
((S^)') = Y.sPis)s{s + 1), we find ((S,)') = 1.62 and 
((Sj)2)o = 1.41 for I^-qP*) and I^-q), respectively, consid- 
erably smaller than ((8^)^)^=3/2 = 15/4 for s = 3/2 and 
{{Sif)s=i = 2 for s - 1. 

Now, we turn to the SDW. The local charge and spin 
distributions shown in Fig. |3| change only slightly from 
the PM to the SDW phase, as also seen in DMFT 
The exchange splitting. A, is within 10% the same for 
all five d-orbitals and amounts to A = 0.37 eV for 
(t/, J) = (8,0.6)eV. The dyz, dxy, dszS-r^, and dj.2_y2 
orbitals contribute almost equally to the magnetization, 
which is TO = 0.74/.tB for ([/, J) = (8, 0.6) eV, whereas 
the dxz orbital, whose band hardly disperses in the anti- 
ferromagnetic direction, is twice as polarized as the oth- 
ers. Its exchange splitting is 10% larger. 

On the right-hand side of Fig. [T] we compare the quasi- 
particle bandstructures of the PM and SDW phases for 
iU, J) = (8, 0.6) eV. Where the PM_E(k)-bands cross the 
E{\i + q) bands (q = [F Y] = [XM]) they split approx- 
imately by A times the overlap of their d-orbital char- 
acters This leads to a noticeable rearrangement of 
the bands. In GT also band-shifts are permitted. Never- 
theless, the GT quasi-particle bands and FS in the SDW 
phase are similar to the results of Stoner theory using 
the PM DFT bands from Hq in ((U and the exchange- 
interaction constant /stoncr adjusted to give the same 
moment 81. 



FIG. 4. (Golor online) GT ordered magnetic moment as a 
function of U for various values of J/U. Our characteristic 
point {U, J) = (8, 0.6) eV is a filled square. Full symbols: 
full atomic Hamiltonian; open symbols: restricted atomic 
Hamiltonian with density-density interactions only. Inset: 
Ground-state phase diagram as a function of U and J. In 
the SDW metal, we also show lines for constant magnetiza- 
tion, m = OAfiB and m = 0.5/iB, respectively. 



The dependence of the ordered magnetic moment on 
the atomic parameters {U, J) in GT is shown in Fig. S) 
We find a remarkably broad region in the {U, J) param- 
eter space with a small ordered moment m < 0.8/iB, in 
contrast to a description based on six Fe electrons in five 
orbitals [iSj, [ij]. As for the five-band model [13J, 



the 

PM metal is the ground state for not too large values of 
the Hund's-rule coupling J, in contrast to Hartree-Fock 
theory where symmetry breaking occurs for small U. 

The atomic Hamiltonian is frequently approximated by 
density-density interactions only, e.g. in DMFT [13, Em 
in order to keep the numerical problems under control. 
One of the magnetization curves in Fig. |4| compares the 
GT result of this approximation with the GT for the 
full Hamiltonian. We see that the former approximation 
largely overestimates the magnetization. 

The inset of Fig. |4| shows that the onset of the SDW re- 
quires a finite J > JciU) « /c-0.017f7 with/c = 0.52 eV 
([/ > 3 eV). Moreover, in the SDW region, the moment is 
to a very good approximation a function of merely the lin- 
ear combination J + 0.017J7. Thus, we then can identify 
I{U, J) = J + Q.017U as the effective low-energy scale for 
the magnetic excitations of our eight-band model in GT. 
That not only the onset but the entire GT magnetization 
collapses onto one curve when plotted as a function of / 
is seen in Fig. |S] This shows that the Hund's-rule cou- 
pling and not the Hubbard-repulsion is the controlling 
parameter for the magnetism in this system. This may 
explain why DMFT studies of LaOFeAs obtain a small 
ordered moment using quite different J7- values lol- 12|. 

Using the PM DFT Hamiltonian Hq, eq. we cal- 
culate the magnetization, m(A), as a function of an ex- 
ternal staggered field A. The resulting A/to(A) = Xdft 
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FIG. 5. (Color online) Ordered moment as a function of / = 
J + 0.017L'' for various values of J/U in GT (diamonds), and 
the inverse susceptibility in DFT Xdft ~ A/m(A) (dots). 
The ordered moment given by spin-DFT (GGA) is reproduced 
by Jdft = 0.82 eV 



is a pure bandstructure function, which in Stoner theory 
provides the selfconsistent magnetization as the solution 
of the equation x5ft("^) ~ ^stoner @]- This function is 
plotted in Fig. [5] As seen from the size of the magnetiza- 
tion jump, the transition from the PM to the SDW state 
in GT is driven by the same band mechanism. Details 
differ due to different underlying PM bandstructures. 

The existence of a sharp 'nose' at small moments is 
caused by the nesting of the hole and electron sheets with 
common xz character, while the strong increase of the 
moment, once m > 0.75/iB in GT, occurs when the xy ex- 
change splitting is so large that the spin-f minority band 
empties at F = Y, and thereby starts to contribute to the 
magnetization. At this point also the zz magnetization 
picks up in GT. This is mainly due to zz hybridization of 
the empty XY minority band near Y. Nevertheless, over 
the entire parameter range, nixz > irixy > riizz > rnxY > 
niyz in GT. For DFT the moment increases faster after 



the nose and rrixz = irixy > mzz > rnxY > myz- 



The 



first difference is due to the steeper bands and the latter 
to the slightly different band positions. Roughly speak- 
ing, the GT of correlated quasi-particles supports the 
Stoner picture of band magnetism for the small-moment 
phase in LaOFeAs, provided that the DFT exchange- 
correlation kernel is substituted by a weak, low-energy 
Stoner parameter. 

The correlated band picture can be further scruti- 
nized from the analysis of the GT susceptibility xgt = 
Agt/™(Agt) (not shown). Following Landau Fermi- 
liquid theory, we write it as xgt/xdft = Qd/i'^ + F^) 
where qd is the average mass enhancement factor and Fq 
is a Landau parameter. The enhancement of the suscep- 
tibility is mostly due to qd{U, J), i.e., XGT/ld depends on 
the interaction through /(m) only, thus explaining why 
the three GT curves in Fig. [S] can be rescaled to a uni- 



versal curve even though their mass enhancements are 
very different. For all m we find that XGT/(9dXDFT) = 
1/(1 + F(f (to)) « 0.8. A small positive value Fl}{m) « 0.2 
reflects the tendency to anti-ferromagnetic spin align- 
ment. The applicability of Fermi- liquid theory supports 
the DFT picture of band magnetism in LaOFeAs. 

In this work, we have employed Gutzwiller theory to 
calculate the quasi-particle bands for a two-dimensional 
eight-band Hubbard model for the valence electrons in 
the iron-arsenic planes of LaOFeAs. In a large region of 
the (?7, J) parameter space, we find a spin-density wave 
ground state of quasi-particles with a small magnetic mo- 
ment m < 0.8/iB, as observed experimentally. When 
charge fluctuations are suppressed by the Hubbard-C/, 
the Hund'-rule coupling J is the effective low-energy scale 
which drives the transition to metallic band magnetism. 

In view of the qualitatively correct DFT description of 
the paramagnetic and also the spin-density wave phases 
(allowing for a renormalization of /dft ^ 0.8 eV to 
^Stoner ~ 0.6 eV), the DFT bandstructure provides a rea- 
sonable starting point for investigations of the iron-based 
superconductors. The Gutzwiller theory offers a micro- 
scopic Fermi-liquid description for the correlation correc- 
tions due to the electron-electron interaction. 

This work was supported in part by the Deutsche 
Forschungsgemeinschaft (SPP 1458 and Bo-3536/1). 
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